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The measurement of higher order cumulants of the current noise generated by a nonlinear meso- 
scopic conductor using a Josephson junction as on-chip detector is investigated theoretically. The 
paper addresses the regime where the noise of the mesoscopic conductor initiates activated escape 
of the Josephson detector out of the zero-voltage state, which can be observed as a voltage rise. It 
is shown that the deviations from Johnson-Nyquist noise can mostly be accounted for by an effec- 
tive temperature which depends on the second noise cumulant of the conductor. The deviations 
from Gaussian statistics lead to rather weak effects and essentially only the third cumulant can be 
measured exploiting the dependence of the corrections to the rate of escape from the zero-voltage 
state on the direction of the bias current. These corrections vanish as the bias current approaches 
the critical current. The theory is based on a description of irreversible processes and fluctuations 
in terms of state variables and conjugate forces. This approach, going back to work by Onsager and 
Machlup, is extended to account for non-Gaussian noise, and it is shown that the thermodynamically 
conjugate force to the electric charge plays a role similar to the counting field introduced in more 
recent approaches to describe non-Gaussian noise statistics. The theory allows to obtain asymptoti- 
cally exact results for the rate of escape in the weak noise limit for all values of the damping strength 
of the Josephson detector. Also the feedback of the detector on the noise generating conductor is 
fully taken into account by treating both coupled mesoscopic devices on an equal footing. 
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I. INTRODUCTION 

Traditional noncquilibrium thermodynamics assumes 
Gaussian fluctuations of the gross variables about their 
mean values^ This assumption is a natural consequence 
of the central limit theorem implying small fluctuations 
of additive variables distributed in a Gaussian way. In 
the last decade there have been extensive theoretical 
efforts^ to calculate deviations from Gaussian statistics 
for electronic current fluctuations of mesoscopic devices. 
The complete knowledge of the number of charges trans- 
ferred through the device in a given interval of time is 
referred to as full counting statistics (FCS). It has turned 
out that FCS reveals details on microscopic processes in 
the device that are not available through mere measure- 
ments of the mean current and the noise variance. This 
can already be seen from a simple example known since 
a long time^ The FCS of a tunnel junction is Poissonian 
when the applied voltage is large compared to the tem- 
perature (eV ^> ksT). In this case charges essentially 
only tunnel from source to drain, and the Poissonian 
statistics points to statistically independent transfers of 
discrete charges. 

In contrast to the substantial literature on theo- 
retical predictions for FCS there are only rather few 
experiments£&£&2ii!ii that have measured deviations 
from Gaussian noise. This is a consequence of the fact 
that these deviations are typically small and require so- 
phisticated experimental techniques to be detected. The 
pioneering work by Reulet et a/£, has measured the 
third cumulant of the noise produced by a tunnel junc- 
tion. Since the noise was measured by room tempera- 
ture electronics, the signal had to be transmitted from 



the cryostat to the amplifier by coaxial cables. There- 
fore, in view of impedance matching, this set-up works 
well only for noise generating devices with resistances 
of order 5011. The more recent experiments employ on- 
chip noise detectors, either quantum point contacts^ or 
Josephson junctions!*^ A first suggestion to use Joseph- 
son junctions as threshold noise detectors was made by 
Tobiska and Nazaro\^ in 2004, and since then vari- 
ous aspects of this idea have been analyzed by several 
authorsffi^i^^^^ 

Two recent experiments§*i2> have studied the noise gen- 
erated by a tunnel junction through measurements of the 
switching rate of an on-chip Josephson junction out of 
the zero- voltage state. The skewness of the noise can 
be extracted from the asymmetry of the switching rate 
with respect to the direction of the bias current. In the 
region of noise activated escape, relevant for the experi- 
ments, the switching of a Josephson junction noise detec- 
tor has been investigated in two recent papers. The work 
by AnkcrholcU^ describes the dynamics of the Joseph- 
son junction in terms of a Fokker-Planck equation driven 
by external noise. An approximate analytical expression 
for the switching rate is obtained for the entire range of 
damping parameters. The subsequent work by Sokho- 
rukov and Jordan^ 7 - employs a path integral formalism 
and accounts for the feedback of the noise detector on 
the noise generating device. The authors also derive 
asymptotically exact results for the switching rate in the 
weak noise limit, however, only for the cases of vanishing 
damping and strong overdamping. In these limiting cases 
the problem simplifies considerably, since the number of 
relevant state variables is halved. The experimentally 
significant parameter range is at intermediate damping. 
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The aim of the present work is to provide, for the 
region of activated escape in the weak noise limit, an 
asymptotically exact solution for the switching rate of 
a Josephson junction in presence of a device that gener- 
ates non-Gaussian noise. The mutual influence of the two 
mesoscopic devices, Josephson noise detector and noise 
generator, will fully be taken into account by treating 
them on an equal footing. Furthermore, the entire range 
a damping parameters of the Josephson junction will be 
covered. 

The article is organized as follows. Sec. [TT| briefly re- 
views a simplified version, sufficient to the present pur- 
poses, of the path integral representation of nonequi- 
librium thermodynamics in terms of thermodynamically 
conjugate variables. This approach was introduced more 
than fifty years ago by Onsager and MachlupiS for the 
linear range near equilibrium and was then extended to 
the nonlinear range by Grabert, Graham, and Greeni 20 i 21 
The method, which is based on the conventional concept 
of Gaussian fluctuations, will then be applied in Sec. ITTTl 
to the thermal escape of a Josephson junction driven 
by Johnson-Nyquist noise. These two introductory sec- 
tions will also serve to introduce the relevant notation. 
The model described in Sec. IIIII will then be extended 
in Sec. IIVI to account for non-Gaussian noise generated 
by a nonlinear device. Finally, Sec. [V] discusses concrete 
results for the experimentally relevant range of parame- 
ters and presents our conclusions. Some more technical 
details are moved to appendices. 



II. PATH INTEGRAL REPRESENTATION OF 
FLUCTUATIONS IN NONLINEAR 
IRREVERSIBLE PROCESSES 

Einstein 2 ^ and Onsager— have related the stochas- 
tic theory of spontaneous fluctuations about equilibrium 
with the deterministic theory of irreversible processes. 
Perhaps the most seminal expression of this relation be- 
tween irreversible processes and fluctuations is the path 
integral representation for the transition probability be- 
tween two macroscopic states. This functional, which 
gives a generalization of the Boltzmann probability dis- 
tribution to the time domain, was introduced by Onsager 
and MachlupA- for the linear range near equilibrium and 
extended to nonlinear processes by Grabert, Graham, 
and Green i 20 ' 21 

Originally, the theory was formulated for closed sys- 
tems where the entropy is the appropriate thermody- 
namic potential. Here we want to apply the method to 
describe mesoscopic systems that exchange energy with 
a cryostat. The modifications are, of course, well-known. 
The entire closed mega-system is divided into the system 
of interest and the heat bath at constant temperature T, 
and the Helmholtz free energy F becomes the relevant 
thermodynamic potential to characterize the system of 
interest. When the state of this system is described in 
terms of the state variables a = (ai, . . . , q,n), the Onsager 



transport equations take the form 

a\ = fi = ^ iijAj , 



(1) 



where the Lu are the Onsager transport coefficients, 
while the 



IdF 



(2) 



are the thermodynamic forces. The transport equations 
are nonlinear, if the thermodynamic forces arc nonlinear 
functions of the state variables a or if the transport coef- 
ficients Lu depend on the state variables. As will be seen 
below, for the problem addressed here, the state depen- 
dence of the transport coefficients is not relevant, and it 
will therefore be assumed that the Lu are constant, they 
may depend on temperature and other external param- 
eters th ough. This simplifies the general theory treated 
in Refs. [20II2TI quite considerably^ 

The state variables a\ can be chosen to be either even 
or odd under time reversal 



f , for even variables 



a\ = e\a\ , sj 



(3) 



— f, for odd variables . 

The Helmholtz free energy is an even variable 

F(a) = F(a) , (4) 

and the transport coefficients obey the reciprocal rela- 
tions 

L L j(a) = eiejLu(a) . (5) 
The matrix Lu may be split into a symmetric part 

(6) 



1 



Du = - [L u + Lji] 



and an antisymmetric part 



1 



Au — ^ [-ku — An] 



(7) 



This implies a decomposition of the deterministic fluxes 
fi into a reversible drift 



n 



(8) 



with the symmetry n(a) = — ein(a), and an irreversible 
drift 



(9) 



with d\(a) — e\di(a). Only the irreversible drift con- 
tributes to the time rate of change of the free energy 



F=-TJ2 Aidi - -T J2 Ai* = -T £ AjAjAj . 



1, .1 



(10) 
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Often, and in particular for the systems treated below, 
some of the state variables do not couple directly to mi- 
croscopic degrees of freedom, and their fluxes are purely 
reversible. We then chose the set of state variables a 
so that the first n variables (oi, . . . , a a , . . . , a n ) are those 
with purely reversible fluxes 



fa 



(11) 



These variables will be distinguished by Greek indices 
a, j3, while the remaining variables (a n +i, . . . , a u . . . , <zn) 
with partly irreversible fluxes will be marked by small 
roman indices i,j. As previously, large roman indices I, J 
run through the complete set from 1 to N. Since the first 
n transport equations of the set |T]) take the form (fTTj) . 
the symmetric parts of some of the transport coefficients 
vanish 



D, 



a,0 



D r 



A,, 



0. 



(12) 



In the stochastic theory of irreversible processes the 
irreversible drift is intimately connected with sponta- 
neous fluctuations about the deterministic motion^ 
These fluctuations can be accounted for by random con- 
tributions ?/i to the thermodynamic forces Ai. Follow- 
ing the approach by Grabert, Graham, and Green ) 20 i 21 
the stochastic theory can be described in terms of a 
Hamiltonian 2 ^ 



H{a, J]) = - ^ DuViVJ + /iW'f ' 
i, J i 



(13) 



which implies equations of motion of canonical form 

dH 
drji 



«i 



/i + ^Aj% 



dH ^ dfj 
oai ^— ' <~> n * 



(14) 



Note that the deterministic transport equations |T]) are 
special solutions of (TT4")) with rji = 0. 

The canonical equations can be interpreted as Euler- 
Lagrange equations and constraints (for the purely re- 
versible fluxes) of an action principle. The action de- 
termines the probability of a fluctuation path, and the 
transition probability from an initial state a(0) = a to a 
final state a(t) = a' may be written as a path integral 



pt(a'\a) = j D[a,rj\ exp |-^-vl[a,rj]| 



with the action functional 



A[a,rj\= / ds^rjiai - H{a,rj) . 
Jo j 



(15) 



(16) 



Since in view of Eq. ([12]) the Hamiltonian (113)) has 
quadratic terms for the r/i only, the action functional is 
linear in the r\ a which act as Lagrange parameters enforc- 
ing the constraints (jlip . The rji, on the other hand, are 




FIG. 1: Circuit diagram of a Josephson junction with critical 
current I c and capacitance C biased by a voltage source V 
via a resistor R. 



random forces describing fluctuations away from the de- 
terministic motion. The Hamiltonian is quadratic in the 
rji because of the underlying assumption of Gaussian fluc- 
tuations. For mesoscopic systems this assumption may 
not be sufficient and an appropriate extension of the ap- 
proach to incorporate non-Gaussian noise will be given 
in Sec. HVl 



III. THERMAL ESCAPE OF A JOSEPHSON 
JUNCTION FROM THE ZERO- VOLTAGE STATE 

In this section the thermally activated escape of a 
Josephson junction form the zero-voltage stated is re- 
viewed utilizing the approach outlined in the previous 
section. 



A. Transport Equations of a Biased Josephson 
Junction 

The state variables of the Josephson junction are the 
charge Q on the junction capacitance C and the phase 
difference (p between the order parameters of the super- 
conductors on either side of the tunnel barrier. The 
time rate of change of the phase is related to the volt- 
age Vj across the Josephson junction by the Josephson 
relatior 



,27 



(17) 



When a voltage V is applied to a Josephson junction in 
series with an Ohmic resistor R, as depicted in the circuit 
diagram Fig. [T] the electrical current / flowing through 
resistor and junction reads 



I=-(V-Vj)=Q + I c sm( l p). 



(18) 



where the second equality follows with the help of Joseph- 
son's relation I s = I c sin(yj) for the supercurrent I s across 
the junction. Combining Eqs. flTSJl with Vj = Q/C, 
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we readily find the deterministic equations of motion 

• - ^9. 

Q = | (v-Q\-icM<p)- 

Clearly, ip is a variable with purely reversible flux 
Let us introduce the free energy 

I c cos(<p) 



(19) 



F(Q^) = F (T,F) + g-A 



V 



and the thermodynamic forces 
1 dF 1 h 



X 



Q 



T dip 

1 dF 

TdQ 



T2e 
TC 



sin(<^>) 



V 
7? 



(20) 



(21) 



The equations of motion (|19p can be then written in On- 
sager form 

2eT 



•■P 



o -■- 



h I \ ML T. 

/ \ ft 7? 



T 
R 




(22) 



Following the approach outlined in the previous section, 
and denoting the conjugate variables to (ip, Q) by (/x, A), 
the Hamiltonian of the system is found to read 



H(cp,Q,r),X) 



-I c $m(ip)\ + 



A 



(23) 



leading to the canonical equations 
dH _ 2e Q 
d/7 ~ Tc 
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Q 



X 



dH 
~d\ ~ 

dH 

dip 

dH 



Uv-9. 

R\ C 



I c sin(y;) 



T , 



= I c cos(ip)X 
1 



:A 



2e ju 



(24) 



While the purely reversible flux <p> remains unchanged in 
the stochastic theory, the flux Q is now supplemented by 
a current (T/R)X describing Gaussian Johnson-Nyquist 
noise from the Ohmic resistor. 



B. Decay of the Zero- Voltage State 



It is convenient to introduce the dimensionless bias cur- 
rent 



V 



(26) 



Then, for < s < 1, the potential has extrema in the 
phase interval [0, 27r] at 



Vweil,to P = arcsin(s) = — =F 8 , 



where for 1 - s « 1 



8 pb V2(l - s) . 



(27) 



(28) 



When the Josephson junction is trapped in the state 
Vwell = f — 8, the average voltage Vj across the junction 
vanishes. However, this zero-voltage state is metastable, 
since the well is only a local minimum of the potential 
([23)) . To escape from the well, the junction needs to be 
thermally activated to the barrier top at tptop = § + 8. 
This process will be observed with large probability, when 
the barrier height is small, which is the case when the 
dimensionless bias current s is close to 1. We shall not 
discuss here escape by macroscopic quantum tunneling j2£ 
which occurs at very low temperatures. 

The decay rate follows from the transition probabil- 
ity from (p wc \\ to ipto P as governed by the path integral 
(|15p . The dominant contribution to the functional inte- 
gral comes from the minimal action path satisfying the 
canonical equations (|24)) . Let us first consider the reverse 
process, the relaxation from the barrier top ytop to the 
well minimum ip we i\. In this case the most probable path 
is the deterministic path, that is a solution of the evolu- 
tion equations (|24[) with ji = A = 0. The two remaining 
equations of motion can be combined to read 

^+A^ + / c sinM^. (29) 
There is a solution^ <p re i ax (t) of (|29|) satisfying 

^relax(-oo) = <£> top , (/? rc l ax (+Oo) = (p well , (30) 

which describes the relaxation from the barrier top to 
the well bottom. Since /x and A vanish, this deterministic 
trajectory has vanishing action (|16|) . 

The minimal action trajectory for thermally activated 
escape from the zero-voltage state ip we ii is a solution 
ip esc (t) of the canonical equations (|24p with 

^esc(-oo) = Vwell) (/? sc(+Oo) = ip top . (31) 



As is apparent from Eq. (|20p . the Josephson junction 
moves in the effective "tilted washboard" potential 



v 

I c cos(<p) + —p) 



(25) 



The first two of the canonical equations (|2"4")l combine to 
give 

2e Cv+ 2eR V + I ^ ] = ^-- (32) 
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Now, the ansatz ip csc (t) = fielax(~t) satisfies the bound- 
ary conditions (|3Tj) and also the evolution equation (|32|) 
provided 

Aosc(i) = — ^V'relaxC— *) = ^<A5Sc(*) , (33) 

where we have used the fact that y re iax(^) is a solu- 
tion of Eq. (|29|) with boundary conditions ([30]) , and that 

^osc(t) = -<£relax(-*), <Pesc(*) = &elax(-*)- The last 

equation of the set (|24|) then gives 



ftC 
~2lT 

ft 

eT 
ft 

eT 



A,, 



(*) 



ft 
2e 



ft 1 



c(-i) 



- 7 c sin (WelaxC-*)) 



(34) 



where we have again used the equation of motion (|29p sat 
isfied by (p Te iax{t) to derive the last line. Now, Eqs. l|33l 
and (|34[) combine to give 



/^csc (t) = 7 C COS ((fesc{t)) X esc {t) 



(35) 



so that the remaining equation of the canonical set of 
equations (|24|) is also satisfied, and the ansatz ip CS c(t) — 
yreiax( — t) gives indeed the minimal action escape path. 

To determine the action (JTHJ) of the escape path, we 
first note that the Hamiltonian (f2"3")l . which is conserved 
along a solution of the canonical equations, vanishes on 
the escape path, since A esc (±oo) = /i csc (±oo) = 0, as can 
be inferred from Eqs. (03]) and [jgl]). Thus 



A — 



dt 



•^esc (^) Qcsc \t) ~t~ /^csc 



ft 

eT 



1/ 

— - 7 c sin(i/3 esc (i)) 



(t) , (36) 



where we have used the first of the canonical equations 
(gill as well as Eqs. (J33J and ^ to express Q csc {t), 
X esc {t), and Hesc(t) in terms of <y9 csc (i). The result 
may now be transformed to read 



A., 




= y [^(Vtop) - E%well) 



(37) 



where the last expression follows from the boundary con- 
ditions (|3Tj) obeyed by (f csc (t) for £ — > ±oo. 

The rate of escape T from the metastable well may be 
written as 



(38) 



where the exponential factor _B is determined by the ac- 
tion of the most probable escape path <^ csc of the path 
integral; 30 Introducing the barrier height 



AU = U(<p top ) - U(<p we a) , 



(39) 



we obtain from Eqs. (|15|) and (|37j) for the exponential 
factor 



B = 



AU 



(40) 



which is just the standard Arrhenius factor for thermally 
activated decay. The pre-exponential factor / requires 
an analysis of the fluctuations about the minimal action 
path and will not be addressed here. 



IV. JOSEPHSON JUNCTION DRIVEN BY 
NON-GAUSSIAN NOISE 

So far we have studied a biased Josephson junction 
driven by Gaussian thermal noise. We now address the 
question how the rate of escape T from the zero-voltage 
state is modified by the presence of non-Gaussian noise. 
To be specific, we shall consider the shot noise generated 
by a normal state tunnel junction, since this case has 
been examined in recent experiments- 9 ' 10 However, the 
theory likewise applies to other noise generating devices 
with short noise correlation times. 



A. Hamiltonian for Non-Gaussian Noise 

Let us consider a Josephson junction with capacitance 
C and critical current I c driven by two noise sources, see 
Fig. [21 A bias voltage Vb is applied to one branch with an 
Ohmic resistor Rb in series with the junction. This part 
of the set-up corresponds to the model treated in the pre- 
vious section. A second voltage Vn is applied to another 
branch with a tunnel junction of resistance Rn again in 
series with the Josephson junction. Experimental set-ups 
are typically more sophisticated, but the circuit diagram 
in Fig. [2] captures the essentials of a Josephson junction 
on-chip noise detector. The current / flowing through 
the Josephson junction is given by 

i = Yr^Il + Ye^Xl = q + Ic sin(ip) . (4i) 



Rb 



R 



N 



Proceeding as in Sec. IIII1 one readily obtains the deter- 
ministic equations of motion 



2eQ 

H~c 

1 



(42) 



Q 



Rb V ^^J + ^r-§)- /cSinM - 

Since the flux (p is purely reversible, the Hamiltonian 
H (ip, Q, r], A) will depend on the conjugate variable rj only 
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FIG. 2: Circuit diagram of a Josephson junction with critical 
current I c and capacitance C biased in a twofold way. The 
branch to the right puts an Ohmic resistor Rb is series with 
the junction and is biased by the voltage Vb- The branch to 
the left is biased by a voltage Vn and Rn is a noise generating 
nonlinear element, specifically a normal state tunnel junction 
with tunnelling resistance Rn- 



linearly, while the dependence on A comprises linear and 
nonlinear terms. In contrast to the case studied in the 
previous section, the nonlinear terms in A will not be 
just quadratic, since the noise generated by the normal 
state tunnel junction is non-Gaussian. As the voltage 
V N = Vn — Vj across the tunnel junctions grows relative 
to ksT/e, the noise generated by the tunnel junction 
crosses over from Gaussian to Poissonian statistics. For 
the current In through the tunnel junction one has 3 - 



IT ) - V » 

{lN) R~n~ 

(SI N (t)6I N (t r )) = C 2 5{t-t') (43) 

(SI N (t)SI N (t')SI N (t")) = C 3 S(t-t')6(t'-t"), 



where Sixty) = In — (In) and 



C 2 = 



^coth 
Rn 

e v N 



2k B T 



Rn 



(44) 



There are of course higher order noise cumulants, but, 
as we shall see, these are not important in the region of 
noise activated switching of the Josephson noise detector. 
The skewness of the noise described by C3 leads to a 
cubic term in A. Neglecting terms of fourth order, the 
Hamiltonian takes the form 31 



H{ip,Q,tj,X) -- 

e(y N -%) 



T 



1 



_y2 _| ^_ ly B _Q_\y y 

2Rb Rb V C ■ 



Ak B R N 



coth 



fcv - %) 
2k B T 



-Rn~[ Vn ~C 



-I c sin(yj)A 



2eQ 



o(x 4 



(45) 



An expansion of the Hamiltonian in powers of A is jus- 
tified, provided the dimensionless quantity e\/k B *C 1. 
As discussed in detail in App. [5] the size of the random 
forces A causing the escape is proportional to the size of 
the fluctuations of the voltage Vj across the Josephson 
junction, and e\/k B is in fact very small, if the decay 
of the zero- voltage state occurs in the region of noise ac- 
tivated escape. Since Vj = Q/C and A are effectively 
proportional to each other, it does not make sense to 
keep higher order terms in Q/C, rather, the two small 
parameters, e\/k B and Q/CVn, should be treated on 
an equal footing. Hence, the term in the second line of 
Eq. (|45|) . which is already of second order in A, can be 
expanded to first order in Q/CVn- Likewise the Q/CVn 
dependence of the term of order A 3 can be dropped. We 
then find 

H(ip, Q, V , A) = H 3 (<p, Q, V, A) + H 3 (<f, Q, 77, A) , (46) 
where 



H2((fi, Q,r), A) 



f J_ (\n_ \ .2 
\2R B 4fc B " 



bias 



R \\ c 



I c sm{tp)A + Yc fJ " 



(47) 



describes Gaussian noise. Here we have introduced the 
bias current^ 

V B , V N 



-'bias — 



Rb Rn ' 



the second noise cumulant 
eV N 



2,N 



R 



coth 



eV, 



N 



2k B T 



and the parallel resistance 
I I 



1 



R 



Rb Rn 



(48) 



(49) 



(50) 



The term 



24k% 



f dC 2 , N Q 
Ak B dV N C 



with the third noise cumulant 



C: 



3,N 



, 2 Vn 
" R 



N 



(51) 



(52) 



includes the leading order effects of non-Gaussian noise. 
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B. Minimal Action Escape Path in the Nearly 
Gaussian Regime 

In the range of parameters studied here, the third or- 
der Hamiltonian ([51) will describe weak corrections to 
the dynamics governed by the Hamiltonian (|47|) . In fact, 
this latter Hamiltonian is precisely of the form of the 
Hamiltonian (f2"3")) studied in Sec. lIIII for a Josephson junc- 
tion in parallel with on Ohmic conductor, provided we 
replace R by the parallel resistance Rn, the current V/R 
by the proper bias current /bias, and T by the effective 
temperature 



7^ — R\ 



R 



T 


C2,N 




Rb 


2k B _ 




T 


eV N 




Rb 


h 2k B R 


N 



coth 



eV N 
2k B T 



(53) 



For eViv *C k B T the tunnel junction generates approx- 
imately Johnson-Nyquist noise and the effective tem- 
perature coincides with the cryostat temperature. On 
the other hand, for eV^ k B T, the tunnel junction 
is a source of shot noise with a noise power propor- 
tional to Vjv- The Josephson junction reacts to the ad- 
ditional Gaussian noise in the same way as to an ele- 
vated temperature.— Approximate expressions for T c ff 
have been presented previously 16 ^ 7 Experimentally, T e g 
can be substantially larger than T. 

The rate of escape T from the zero- voltage state of the 
Josephson junction will again be of the form (I38p . where 
the exponent B now takes the form 



B = Bn 



(54) 



where we have made use of Eq. (I53|) , and 
fi = --Q— = I c cos(tp)\ 



1 2e fi 

A — 



dH 

dQ R\\C" hC 

1 dC 2 , N \ 2 
Ak B dV N C 



(57) 



Now, the two differential equations (|56[) of first order can 
be combined to one second order differential equation 

ft 1 T 
—Op + — — •p + Ic sin(ip) = /bias + -^-A + h , (58) 
le Ze /t|| K\\ 



where 



~2~ C3,n^ 2 



ft dC 2iN . . 
4e/c B dVj 9 



N 



(59) 



is the additional noise current arising from i/3 . Likewise, 
the Eqs. |57|) combine to give 

^C\- ^^-A + / c cos(^)A = I 3 \, (60) 



2c 



where 



2e R v 



j, _ ft 8C 2 ,n ; 
3 ~ 4ek B dV N 



(61) 



again results from H 3 . 
We now make the ansatz 



with 



<p esc (t) = p 2 (t) + tp 3 (t) 



B 2 = 



AU 

k B T c g 



(55) 



The exponential factor B 2 is determined by the action of 
the approximate escape path (p 2 (t) that solves the canon- 
ical equations of motion resulting from the second order 
Hamiltonian (|4T)l . The second cumulant (|4"9"|) of the noise 
generated by the normal state tunnel junction is taken 
into account in terms of the effective temperature T e ff. 
To include the effects of the third cumulant 63, at, one 
needs to determine the deviation ip3 (t) of the escape path 
from <p 2 (t). To this purpose we start with the canonical 
equations that follow from Eqs. (g7]) and d3TJ) . We 

find 



2eQ 



dH 
dH 



+ 8kJ Cs ' NX 



R\\ C 



- I c sin(</j) 
1 dC 2 , N Q 



/off 

~Rn' 



2k B dV N C 



A, 



(56) 



(62) 



Acsc(i) - A 2 (i) + A 3 (i), 



where <p 2 (t) and X 2 (t) are the solutions of (|58[) and 
for / 3 = Jg = 0, while <£z(t) and A 3 (i) describe the mod- 
ifications of the path arising for finite I3 and I 3 . For 
I3 = 0, the equation of motion (|55|) is of the form of the 
evolution equation (j32|) studied in Sec. IIII1 and we can 
proceed as there. Provided s = /bias//c < 1, the po- 
tential U(ip) = —(h/2e)I c [cos((p) + sip] has a minimum 
</? wc ii and a maximum <pt op in the phase interval [0, 2ir]. 
From a solution <Preiax(i) satisfying 



ft_ 

Ye 



Cip- 



2~e~R 



(p + I c sin(^) 



' bias 1 



(63) 



and the boundary conditions (|30[) . we obtain an escape 
path satisfying the Eqs. (JSHJ and ([HO]) for / 3 = I 3 = 
and the boundary conditions 



tp 2 (-oc) = tp wc n , tp 2 (+00) = tp top 



A 2 (-oo) = , A 2 (+oo) =0 



(64) 
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by putting 



lf 2 {t) = <y2rclax(-£) 

h 



eT cS 



02 (t). 



(65) 



Next, we insert the ansatz (1621) into the evolution equa- 
tions (|58f and (|60|) and keep only terms that are linear 
in the quantities <p3, A3, 23, and 1% which describe cor- 
rections to the Gaussian case. Taking advantage of the 
equations of motion satisfied by Lp 2 and A2, we obtain 



2e 
and 



h 1 
Ye~R\ 



(p 3 + I c COs((/? 2 )c/5 3 



111 

R 



II 



A3 + /3, (66) 



— CA3 — — A3 



2e 



2eRu 



J C C0S(</3 2 )A3 



err 



-ip 2 [7 c sin(^ 2 )^3 + 



(67) 



where I3 and I' 3 defined in (l59|) and (fBTj) are now evalu- 
ated with the leading order solutions (|65|) . Hence 



and 



2(fc s T eff ) 2 V2e 



C: 



3,N 



2k B T cS - 



dC- 



2,N 



dV N 



■ 2 

^2 ■ 



knT, 



2e 



<9G 



2,iV 



JV 



■ <^2 ■ 



(68) 



(69) 



We shall see that an explicit solution of these evolution 
equations is not required to determine the action. 



C. Action of Escape Path 

Since the Hamiltonian (|4"6")) vanishes along the escape 
path, the action may be written 



dt 



•^cscQcsc AescV-'esc 



(70) 



where we have made a partial integration with respect 
to the first line of Eq. (|3"6")) . From Eq. (fSo) , we have 
Qcsc = {h/2e)C(pesc, while Eq. (|57|) implies fi csc = 
I c cos(</?esc)A eS c- Inserting this as well as the ansatz (|62[) 
into the action ([70]) . we find after disregarding terms of 
second order in ip^ and A3 



where 



Ao = / fit 



A csc = A 2 + A 3 , (71) 



■^-C\2(f>2 - i c cos(tp 2 )A2^2 , (72) 
ze 



and 

A3 



dt 



2e 



C (A 2 (/3 3 + A 3 c/3 2 ) 



(73) 



-7 C cos(tp 2 ) (A 2 <^3 + A3CP2) + sin(( ) C2)A2^ 2 <P3 



Now, the deviations (f3 and A3 from the path of the Gaus- 
sian model are caused by the currents I3 and Jg given in 
Eqs. ((68)) and lj69|) . These currents depend on the third 
noise cumulant C^jv and on the derivative BCz^n / 9Vn 
of the second cumulant. The detailed evaluation of the 
action in App. [Bl shows, that these two factors influence 
the action A3 only in the combination 



C3 = Q 



3,N 



3k B T, 



dC- 



2,N 



cff" 



dV, 



N 



(74) 



A corresponding reduction of the effective third cumulant 
was already noted by Sokhurokov and Jordan^ for the 
limiting cases of weak and strong damping. The second 
term in Eq. I|74[) arises from the feedback of the Joseph- 
son junction on the noise generating junction, which is a 
consequence of the finite voltage Vj that builds up dur- 
ing escape. Experiments are usually done in the regime 
eVjv S> k B T, where 



C 3 



C: 



3.N 



1 



3k B T c s 
eV N 



C: 



3,N 



1 - - 



3Rb 

2~R b +Rn 



2k B T 
eV N 



R 



N 



(75) 



so that the feedback becomes negligible for Rn ^> Rb- 
In the opposite limit the feedback even changes the sign 
ofC 3 . 

As shown in App. [Bj repeated use of the equations 
of motion satisfied by <p 2 , A2, ^3, and A3 allows one 
to express A3 entirely in terms of f2(t)- By virtue of 
Eq. I|65[) . tf2{t) is time reversed to the deterministic tra- 
jectory i^relax(i) describing the relaxation from the bar- 
rier top. Accordingly, the result (|B20[) in App. [Bl may be 
written as 



A, 



2k f 



{k B T eS y V2e 



C 3 J- 



where 



1 r°° 

J =-qJ 



(76) 



(77) 



Thus, the non-Gaussian correction to the rate exponent 
([541 reads 



Bo 



1 



{k B T cS ) 



C 3 J- 



(78) 



What remains to be determined is the quantity J, which 
describes a property of the system in the absence of noise. 
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Let us introduce the energy function 



DISCUSSION 



E M) = 1(J£) CV 2 + U{^), (79) 

where U(<p) is the potential ([2"5)) with V/R replaced by 
-^bias = sl c . The time rate of change of E reads 

d f K \ ^ h 

Jt E= \Ye) C(p( ? + 2~e [Ic Sin( ^ " /bias] ' (80) 



which, using the equation of motion (|63p satisfied by 
( / 3 reiax(^)j may be written as 



dt 



nY l 



2eJ R 



!! 



(81) 



Along the deterministic trajectory ip tc \ ayi (i) we may look 
upon E as a function of </?. Then 



dE _ idE _ ay 1 . 

dt \2e/ ' 



and from Eq. (|T9"|) we have 



(82) 



(83) 



In this section we will give some concrete results in the 
experimentally relevant range of parameters. 



A. Dimensionless Quantities 

It is convenient to formulate the theory in terms of 
dimensionless quantities. Introducing the plasma fre- 
quency of the Josephson junction at vanishing bias cur- 
rent 



2el c 



p ~ V n c 

the result ([So]) may be written as 

J = 3 , 

where 



(87) 



1 rfwell 

j = -- / d<p(e-u) 

6 JVtap 



(89) 



which combines with Eq. ([82)1 to yield 



dE h 1 2 , 



dtp 



2eR\\ V C 



U), 



(84) 



where the sign is determined by the fact that E decreases 
along the trajectory. 

The function E(ip) can easily be determined by numer- 
ical integration of Eq. ([84)1 . One starts from tp = tp t0 p 
with energy E(ip top ) = U(ip top ) and integrates towards 
smaller ip with the + sign of Eq. ([54")) until the first turn- 
ing point with E((p) = U(<p) is reached. There, the inte- 
gration continues towards larger values of ip with the — 
sign of Eq. ([51)1 up to the second turning point and so 
on, until the trajectory ends at E(ip wc n) = U((p we n)- 

By virtue of Eq. ([8"3")l the formula ([77| may be written 

as 



3C U 



dttp(E-U) 



(85) 



Changing from an integration over time to one over 
phase, we get 



J 



1 2e 



3C V h 



dip (E-U), (86) 



Vtop 



where the integration starts at ptop and goes back and 
forth between the turning points until it ends in ip we n- 
The determination of the effect of non-Gaussian noise on 
the rate of escape is thus reduced to an integration of the 
first order differential equation ([84")) and the evaluation 
of the integral 



is a dimensionless integral given in terms of the dimen- 
sionless energy 



2eE 1 



h I c 2lu; 



and the dimensionless potential 



:<p~ + U 



2eU 

u = — = — cos(p) — sip . 

n Ic. 



(90) 



(91) 



From Eq. ([75)1 . the correction to the exponential fac- 
tor of the rate may then be written as 



W (k B T cS f 



(92) 



To determine j from Eq. ([89)) . one needs to solve the 
dimensionless form of Eq. ([84")) , which reads 



de 
dip 



± 7N /2(e- w) 



where 



7 = 



(93) 



(94) 



is the dimensionless damping coefficient, which coincides 
with the inverse quality factor Q = R^Cw p at vanishing 
bias current. 
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B. Strong Damping 

Let us first discuss the limit of strong damping 7> 1. 
The Josephson junction noise detector cannot operate in 
this limit, because after escape from the metastable well 
the phase will be retrapped in the adjacent well of the 
tilted washboard potential, so that only a short voltage 
pulse builds up. Nevertheless, the behavior in this limit 
is instructive, since explicit analytical results can be ob- 
tained. To solve Eq. ([93]) . we make the ansatz 



e = u + k 



and find 



This gives 



dn du 

dip dip 



± 7V2k . 



_ 1 f du dn , 

H = ±—=- -3- + -T- » 

V27 \dip dip, 



(95) 



(96) 



(97) 



so that the dimensionless kinetic energy n is of order I/7 2 
for large 7. The leading order solution 



1 f du 



27 2 \dip 



(98) 



satisfies the boundary conditions e = u, i.e. k = 0, for 
<P = '/'top and ip = ip W eii- Inserting Eq. (|9"8"|) into Eq. 
we obtain 



tp{t) oscillates back and forth in the potential well and 
looses energy only very gradually. Let us consider a seg- 
ment of the trajectory starting at a turning point ip+ on 
the barrier side of the potential, oscillating through the 
potential well to a turning point ip_ on the opposite side, 
and traversing the potential well again to a turning point 
ip' + . From Eq. (|93|) we find for the energy along this path 
segment 



e(ip) = e(ip + ) + 7 / dtpy/2 (e - u) 
J<p+ 



± 7 / (kpy/2 (e-u), (102) 
J ip _ 

where the + sign holds for the oscillation form ip + to 
if- , and the — sign on the way back from ip- to ip' + . For 
7 <C 1, this gives 

r<p- 

e(ip) = e+ + 7 / d(fy/2 (e+ - u) 



± 7 / dip y /2{e+- u)+0( 1 2 ), (103) 
</ 

where e+ = e(y>+) = u{ip+). This result can now be 
inserted into Eq. (f89|) . to find for a segment of the ip- 
integral form ip+ over if- to </>+ 



A; = -1 I | V (e - u) + ^ + d</> (e - «) 1 (104) 



i_ r wci1 / 

6 7 2 L tOD W 



(99) 



In the overdamped limit, there are no turning points, but 
the phase gradually slides down from iptop to (/Veil • Using 
Eqs. (HI]) and l[9T]). Eq. ^ is readily evaluated with the 
result 



J = 



(l + 2s 2 ) arccos(s) - 3s VI - s 2 



6 7 2 



(100) 



Now, the observed escape events occur typically for val- 
ues of the bias current /bias close to the critical current 
I c . Then, l-s« 1 and Eq. (|100p can be expanded to 
yield 



3 



8^ (1 _ s)5/2 ! 



45 



(101) 



7 



This latter formula is in accordance with the result by 
Sukhorukov and Jordan 17 in this limit. 



C. Very Weak Damping 



= -7/ dip dtp'yfi [e+ - u(ip')} + 0( 7 2 ) , 

where we have taken into account that the difference be- 
tween ip + and ip' + is of order 7. 

On the other hand, Eq. (|103[) gives for the change Ae 
of the energy during one oscillation period 

Ae = -2 7 [ V+ dip^/2 [e+ ~ u{ip)} + C( 7 2 ) . (105) 



(106) 



Eqs. (|104p and (|105|1 combine to yield 
^ = -/(e) + 0(7), 

where 



= l /;_ + ^(^ + -^) y^-nfo,) 
3 PJ+ dipy^~ujpj 

Dividing the integral (|89p into segments of the form 
(|104| . we can transform the integral over ip into an inte- 
gral over energy. Using Eq. (|106|) . we then obtain 



Next we consider the case of a very weakly damped 
Josephson junction, i.e., 7 < 1. Then the trajectory 



i(Vtop) 

./ = / def(e) . 



(108) 
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Let us again study specifically the experimentally im- 
portant range 1 — s <C 1. Then, the relevant range of ip 
values lies in the vicinity of \. Putting 



we find for the potential (|9T 



where 



= -l s + V2(l-s) 



? = i> - 



3/2 



(109) 



(110) 



(111) 




With the scaled dimcnsionless energy 

e = -| S + V / 2(l- S ) 3/2 e, 



(112) 



the result (| 1 08[) with (|107p can be transformed to read 



FIG. 3: The scaled dimensionless energy e is shown as a func- 
tion of ip for 7 = 0.25. The energy decreases as the trajectory 
moves back and forth in the potential depicted as a grey 
line. 



de 1 — — — — — ^ , (113) 



where "0- an d "0+ are the negative and smallest positive 
roots of s(ip) = ip — ^ip 3 = e, respectively. The remaining 
integral is just a numerical factor independent of s, and 
a numerical evaluation gives 



j = a(l - s) 2 , with a = 0.79. 



(114) 



This result is in accordance with the findings by Sukho- 
rukov and Jordan— in the limit of vanishing damping. 



proper sign back and forth between the turning points 
until the integration ends at e(— 1) = <r(— 1) = — §. A 
typical solution is depicted in Fig. [31 In scaled units 
Eq. (gU) takes the form 

j = -!{!- a) 2 J (H7) 

where the integral follows the ^-path back and forth be- 
tween the turning points. Since the differential equation 
f| 1 15(1 depends on s and 7 only in the combination 7, we 
put 



j = _(i _ s) 2 Wtf) : 



(118) 



where 



D. Intermediate Damping 

In experiments typical values of the dimensionless 
damping coefficient 7 are small but nonvanishing. The 
factor j in formula (j9"2"|) for B3 must then be determined 
from Eq. (|89"|) using the solution of the differential equa- 
tion ([93]) . While a numerical evaluation is straightfor- 
ward for arbitrary values of s, we shall focus on the exper- 
imentally relevant range 1 — s <C 1. In terms of the scaled 
quantities introduced in Eqs. (fTU9]) - fTH)) , Eq. (S3]) reads 



where 



2 V 



1 - s 



(115) 



(116) 



This differential equation has to be solved with initial 
condition e(l) = ^(1) = |, and integrated with the 



w(7) = -y #0-0 



(119) 



The function W(j) determines the correction B3 of the 
exponential factor of the rate for arbitrary damping 
strength in the range of bias currents close to the critical 
current. 

From Eq. (|114[) . we obtain 



W(0) = 1.188. 



while Eq. (jlOlj) gives for 7 > 1 



W(j) « - - 
15 7^ 



(120) 



(121) 



where we have made use of Eq. (|116j) . In between these 
limiting results, the function needs to be determined nu- 
merically. A list of data points is provided in Table[H and 
the function W(7) is depicted in Fig. 2] together with the 
findings of previous worksj 16 ' 17 This should facilitate the 
comparison with experimental results. 
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7 





0.025 


0.05 


0.075 


0.1 


0.125 


0.15 


0.175 


W 


1.188 


1.185 


1.179 


1.169 


1.157 


1.142 


1.125 


1.107 




7 


0.2 


0.225 


0.25 


0.5 


0.75 


1.0 


1.5 


2.0 


W 


1.087 


1.066 


1.043 


0.797 


0.574 


0.409 


0.218 


0.129 



TABLE I: Some numerical values for W as a function of 7. 



T 




II , 1 , 1 , 1 — 

0.5 1 1.5 

Y 



FIG. 4: W is depicted as a function of 7 (straight line). 
Also shown are the results of Ref. |l7J for vanishing damping, 
Eq. (|120|l . (dot), and in the strong damping limit, Eq. (|121|l . 
(dotted line) . The approximate result of Ref. is depicted 
as a dashed line. 



E. Conclusions 

We have presented a theory for a Josephson junc- 
tion detecting non-Gaussian fluctuations by means of the 
noise driven escape out of the zero-voltage state of the 
junction. It has been assumed that the device is op- 
erated in a regime where the barrier of the washboard 
potential is overcome by activated processes. This is al- 
ways the case if the temperature is not too low and/or 
the junction capacitance is not too small. The study was 
based on the theory of irreversible processes and fluctua- 
tions developed by Onsager and Machlup 19 and Grabert, 
Graham, and Greeni 20 ' 21 An extension of the method to 
account for non-Gaussian fluctuations was outlined^ In 
this approach the random motion of the system is de- 
scribed in terms of the state variables and the conjugate 
forces. The force A conjugate to the electric charge Q, 
which appears naturally in this approach, plays a role 
similar to the counting field introduced in the more re- 
cent approaches to determine the full counting statistics 
of electronic devices 4 

A nonlinear noise generating clement in series with the 
Josephson detector modifies the rate of escape out of the 
zero-voltage state. The main effect comes from the sec- 
ond noise cumulant C2,n- However, this Gaussian part 



of the noise is detected by the Josephson junction in the 
same way as Johnson-Nyquist noise. Therefore, as was 
shown explicitly, the second noise cumulant can be de- 
scribed in terms of an effective temperature T e g. Devi- 
ations from the accordingly modified Arrhcnius law are 
thus due to higher order noise cumulants. The fluctua- 
tions causing the escape from the metastable well lead to 
fluctuations of the voltage Vj across the Josephson junc- 
tion. It has been shown that these voltage fluctuations 
are small compared to fcgTeg/e, which implies that the 
dimcnsionlcss random force eA/fcs causing these fluctua- 
tions is always small compared to 1. Since the n th order 
noise cumulant gives rise to terms of order (eA/fcs)", de- 
viations from the modified Arrhenius law essentially only 
arise from the third noise cumulant C^.n, and these cor- 
rections are typically small. However, the third cumulant 
is odd under time reversal and the sign of the effect de- 
pends on the direction of the bias current. Comparing 
rates for pulses tilting the potential to the right and the 
left, respectively, the correction B 3 can be extracted^^ 
A Josephson junction threshold detector operating in the 
regime of noise activated escape thus can measure the 
third cumulant, the skewness of the noise, only. Another 
effect of the fluctuations of the voltage Vj is a feedback of 
the Josephson detector on the noise generating device as 
described by the effective third noise cumulant C3 defined 
in Eq. (d). 

The modification of the rate exponent due to the skew- 
ness of the noise has been determined for arbitrary damp- 
ing strength of the Josephson junction detector. Thereby, 
the theory developed goes considerably beyond the re- 
sults of previous works ) 16 i 17 that were restricted to lim- 
iting values of the damping strength or based on approx- 
imations. Explicit results where given for the case when 
the bias current is close to the critical current, which im- 
plies that the relevant part of the washboard potential 
can be described by a cubic potential. The effect of the 
skewness of the noise on the rate is, however, larger for 
smaller values of the bias current. Experimentally, the 
range of relevant bias currents can be influenced by the 
form of the applied current pulses. The theory presented 
here can readily also be evaluated for the exact form of 
the washboard potential allowing for results for any value 
of the bias current and all damping strengths. 

To be explicit, we have presented the theory using the 
example of a normal state tunnel junction as noise gen- 
erating device. However, the theory readily also applies 
to other noise generating elements, provided the corre- 
lation time of the noise is much smaller than the period 
of plasma oscillations of the detector. Finally, in this 
article, only the exponential factor of the rate has been 
determined. The corrections due to the skewness of the 
noise were found to be rather small, and they need sophis- 
ticated experimental techniques to be detected reliably. 
Corrections to the pre-exponential factor of the same or- 
der of magnitude are entirely negligible, so that safely 
the prefactor of the standard Gaussian noise theory can 
be employed. 
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APPENDIX A: VALIDITY OF NEARLY 
GAUSSIAN APPROXIMATION 

In this appendix we investigate the range of validity of 
the nearly Gaussian approximation used in Sec. HVl Since 
the leading order term <p 2 (t) of the most probable escape 
path is the time reversed relaxation path ip vc \ ayi {t) , the 
order of magnitude of the phase velocity ip during escape 
coincides with that during relaxation. 

Let us first consider the case of weak damping. The 
trajectory ^roiax(i) starts with vanishing phase veloc- 
ity at the barrier top. The largest kinetic energy 
<fi 2 arises when the potential minimum </? W eii 
is reached for the first time. For weak damping the ki- 
netic energy then almost equals the potential energy dif- 
ference AU. Accordingly, the voltage Vj — (h/2e)(p sat- 
isfies 



Vj < 



2AU 

C 



As damping increases the phase velocity and, accord- 
ingly, the maximal voltage across the Josephson junction 
decreases, so that Vj will never exceed the estimate (|A1[) 
in the entire range of parameters. 

The plasma frequency of the Josephson junction at fi- 
nite bias current 



(A2) 



is the frequency of small undamped oscillations about the 
minimum </? W eii of the potential (|25[) . For i5 < |, which 
is the case for 1 - s < 1, Eqs. - (J25J) yield for the 
barrier height (|39|) of the potential 



AU 



Til, 
~3e 



C X3 



This can be combined with Eq. (|A2[) to give 



Hujp(s) 



2ehI r S 



e< /(»Ar 



(A3) 



(A4) 



C S V c 

The bound (|A1I) for the size of the fluctuations of Vj may 
thus be written as 



eVj < -j= hhj p {s) . 



(A5) 



In the region of thermally activated escape^ one has 
hujp(s) <C fcsT c ff. In view of Eq. (|A5I) this implies 



eVj S 
T off V3 



(A6) 



so that eVj/kgTcfi is a small dimensionless parameter 
along the most probable escape path. 

Now, the leading order contribution A2 to the force A 
causing the escape is determined by Eq. (j6"5)) , entailing 
the estimate 

n . _ 2Vj 

eT e ff T cff ' 
which combines with the inequality (|A6|) to give 



A 



(A7) 



eX 



< 1 



(A8) 



This shows that an expansion of the Hamiltonian in terms 
of A, as done in Eq. (|45|) . is indeed justified. The terms of 
third order in A arc then small, so that ips and A3 describe 
in fact small corrections to if 2 and A2, respectively. 

Because of the weak effects of non-Gaussian statistics, 
the correction B3 to the exponent of the rate is also small. 
From Eqs. (J55]) and Ufi2}, we find 



B3 
B 2 



(2e) 



fcsToff 



AU 



J 



For S -C f we can insert Eqs. (|1 18[) and 
Eq. ([28]). we then find 

2 n 



Tno ri 



B 2 16 \k B T e sJ e 2 I, 



WS. 



(A9) 
Using 

(AlO) 



(Al) 

Hence, the effect of the skewness of the noise vanishes 



proportional to (1 — s) 1 / 2 as the bias current approaches 
I c . The ratio B3/B2 can be seen as a product of three 
factors 

(All) 



B 2 



1 

To 



Hujp(s) 
^sT off 



^3 TI , 

x x W . 



e 2 I c 



where we have made use of Eq. (|A2|) . Now, in the regime 
of activated decay the first factor (1 / '16)(fkv p (s) / 'ksT c s) 2 
is very small, while the last factor W is of order 1 for weak 
to moderate damping. Hence, one needs a large factor 
C3/e 2 J c to get observable effects from the skewness of the 
noise. Since Cz,n is proportional to Vat, this means large 
Vn, in particular, eVjv S> fc^T, so that the estimate (|75|) 
for C3 applies. To minimize the reduction of 63,^ via 
the feedback effects described by Eq. ([75]) . one needs to 
choose a bias resistor Rb well below Bn- Then the factor 

C3 C3,jv Vn 
e 2 L ^ 



e 2 L 



RnI c 



(A12) 



This means that the current Vn / Rn should be large com- 
pared to I c and thus needs to be largely compensated by a 
current Vb I Rb in the opposite direction to keep the junc- 
tion biasing current (|48|) below I c . Experimentally, this 
compensation problem is addressed by employing more 
sophisticated set-upsi*^ 
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APPENDIX B: EVALUATION OF ACTION OF 
ESCAPE PATH 



In this Appendix we evaluate the expressions (|72| and 
(175)) for the action of the escape path in the nearly Gaus- 
sian approximation. Inserting the result (|65|) for A2, one 
obtains from (17211 



2 

T e tt 

2 

T e fi 



dt 



dt 



A) C0202 - -Z-Ic COs(<£ 2 )¥>202 

ze / 2e 



d 1 / h 



dt2\2e 



d h 



(Bl) 



Now, 02 vanishes at the integration boundaries and 
— (h/2e)I c [cos(y>2) + f2 sin((^2)] coincides there with 
U(ipto P ) and U((p wc \\), respecti vely , since sin(y we u) = 
sin((^to P ) = s. Accordingly, Eq. (|B1[) yields 



A, = 



2AU 



off 



(B2) 



which gives the exponential factor (|55|) of the escape rate. 
After expressing A2 in terms of (f2 and putting 



A, = 



eTcff 



-A, 



(B3) 



we obtain from Eq. ([75]) for the leading order non- 
Gaussian part of the action 

A * = ^J dt {l (~) <?(0 2 3 + 2 A 3 ) (B4) 



-i" c [C0S(^ 2 ) (02<^3 + V2A3) - Sin(^ 2 )^2^2^3] 

e 



The integral in the first line gives after partial integration 



.4, 



1 



3, part I 



cir 



dt- 



1 fh 



2 V e 



0203 + <^2A 3 



(B5) 

In this expression we can eliminate the second order 
derivatives 03 and A3 by means of the equations of mo- 
tion (l66|) and ([67]) . Taking the definition (|B3j) into ac- 
count, we get 



.4 



3, part I 



eft J- 



dti (p 2 



h 1 

2e"^ 



-7 C cos(y> 2 )¥>3 + — ^-A 3 + J 3 
ei£|| 



+ <P2 



A 3 - / c cos0 2 )A 3 



+02 (/ c sin((/5 2 )<^3 + I3) 



(B6) 



This result can now be inserted into (|B4|) . After a partial 
integration of the ^2^3 term and a further partial inte- 
gration along the lines I c [sui(ip2) l P2 l P2 — cos(y>2)02]v?3 = 
I c [—(d/dt) cos(<p 2 )<P2]<P3 — > -f c cos(<^ 2 )¥'2<P3 J one obtains 



-4, 



dt \& ( J 3 + (B7) 

e-teff J-00 I 



^"ET^ 2 - 2I C COs(<p 2 )<P2 ) (A3 - 3 ) 

2e K\\ 



From Eqs. ([55]) and (|M|) . we see that 



02 (^3 + Pzl'a) = 



1 / h 



2 V2e7 (fc B T eff )- 



• C 3 ,n 02 (B8) 



1 



fc B T efr \2eJ dV } 



H \ 8C2,N 1 -3 • •• \ 

¥>2 - VWWl) 

/ N 



Now, under the integral <^2 0202 — ^(d/dtjhip^ 
— \(p\, so that 952 (^3 + ^2^3) can be replaced by 



1 / h 



02 (h + V2I3) -»r r 



2 V2e7 (fcsTeff)' 



■C 3 ^2 (B9) 



where 



C3 = C: 



3,iV 



3fc s T c ff 



dC, 



2,N 



dV N 



(BIO) 



Since the action (|B7[) depends on A3 — 03 only, is is nat- 
ural to make the ansatz 



A3 = ^3 + A3 



(Bll) 



From Eq. (|B3[) and the equations of motion (|M|) and (|S7[) 
one then finds 

AcA 3 + ^-Aa 3 + I c cos(^ 2 )A 3 = wis ~ /s ■ (B12) 
ze ze rtii 



Using Eqs. (|68|) and (|69|) . the right hand side may be 
written as 



02/ 3 - I3 = - 



1 



2e/ (fcsTeff)' 



■C 3 2 02, (B13) 



where again the cumulants appear only in the combina- 
tion (|BI0)) . 

We can now employ the evolution equation (|B12|) to 
express the term proportional to I c in the action (|B7[) in 
favor of terms with a purely polynomial dependence on 
(p 2 - Using also Eqs. fB^I . (|B11[> . and (|B13j) . we find 



-4 :i 



cir 



eft 



1 / ft 



1 



2 V2e/ (k B T cl 



72 C 3 



h 1 • »/ „ 
---^3 + 2^2 



AcA 3 + A 1 a 3 

2e 3 2ei? M 3 



1 



C3 0202 



(B14) 
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After partial integrations along the lines tp 2 h/ 3 — > (p 2 h! 3 , 
(p 2 A' 3 -> -0 2 A 3 , and (p 2 (p 2 <f 2 = y>i{d j d£)\<fy -> -5^, 
this simplifies to read 



-4, 



eT c ff 



d« - 



1 



-G</? 2 - 7r"5~V2 



2ei? 



2^ 

a; 



1 



(k B T cS )' 



■c 3 y\ 



(B15) 



Comparing the form of the evolution equation (|B12|) 
with the one satisfied by ip 2 , namely Eq. (|58p for ^3 = 0, 
we are led to the ansatz 



A' 3 (t) = A(t)^ 2 (t) . 



(B16) 



Inserting this into Eq. (|B12|) and using the evolution 
equation for cp 2 as well as Eq. (|B13|) . we find that A(t) 
obeys the differential equation 



2e 



2 A 



R\\C 



A <p 2 



A 



R\\C 



A) fa 



kBT c s 



C3 fa fa 



(B17) 



which is satisfied, provided 



A- 



R\\C 



-A 



1 h 



3C 2e ( kB T eS y 



■C 3 fa. 



(B18) 



When the ansatz (|B16j) is plugged into (|B15j) . we obtain 
a term proportional to A fa fa, which under the integral 
can be replaced by —\A(p\. Accordingly, we find 



M = 



eT G ff 



2\2eJ {k B T oS y 



■c 3 & 



2e 



-CA 



h 1 

Yell 



-A 



(B19) 



Finally, in the integrand, the expression between squared 
brackets can be transformed by means of Eq. (|B18[) to 
yield for A3 the compact result 



A* 



1 



3 \2eJ (k B T cff y 



dt ipl . (B20) 
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